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Static properties of SU (3) multiskyrmions with baryon numbers up to 6 are estimated. The 
calculations are based on the recently suggested generalization of the SU (2) rational map 
ansatze applied to the SU(3) model. Both SU{2) embedded skyrmions and genuine ST/ (3) 
solutions are considered and it is shown that although, at the classical level, the energy of the 
C\ . embeddings is lower, the quantum corrections can alter these conclusions. This correction 

to the energy of the lowest state, depending on the Wess-Zumino term, is presented in the 
most general case. 
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1. Topological soliton models, and the Skyrme model among them have recently 
generated a fair amount of interest because they may be able to describe various properties 
of low energy baryons. So far, most of such studies have involved the SU(2) Skyrme model, 

7— ( ] whose solutions were then embedded into the 577(3) model. This is justified as the solutions 
of the SU(2) model are also solutions of SU(N) models. However, there exist solutions of 
the SU(N) models which are not embeddings of SU{2) fields and it is important to assess 
O \ their contribution. 

As with solutions of the 577(2) Skyrme models, the solutions of the SU(3) model, with 
very few exceptions, can only be determined numerically. Like for the 57/(2) case 0] one 
starts with a harmonic map ansatz which gives fields "close" to the genuine solutions Tp, |[] . 
gn! Then one can use these fields as starting configurations of various numerical minimization 
schemes. For the 577(2) fields such an approach was carried out in 0J where it reproduced 
the results of the earlier numerical approaches []5], [| and so showed that the harmonic map 
approximations are very close to the final fields. The real solutions have energies only few 
hundreds of MeV lower than the harmonic approximations and the baryonic charge and 
energy distributions do not look very different etc. This has justified the use of harmonic 
field approximants in calculating various quantities when estimating quantum corrections 
to the classical results. 

In this paper we take a look at the nonembedding solutions of the SU(3) model. 
First we recall the results of Q and use them in a SU(3) numerical minimization program 
to estimate the values of energies of true solutions. The SU(3) variational minimization 
program has been rearranged for this purpose to allow the consideration of quite general 
field ansatze. This is discussed in the next section. 

In the following section we discuss various quantum corrections and compare our 
results with the similar results for the embeddings. The paper ends with a short section 
presenting our conclusions. 

2. The Lagrangian of the SU(3) Skyrme model, in its well known form, depends on 
the parameters F n and e and is given by || [|: 

C = ^Trl^ + ^Tr[l„ I,} 2 + ^Tr(U + - 2) (1) 



Here U e 577(3) is a unitary matrix describing the chiral (meson) fields, and = U^d^U. In 
the model F n should be fixed at the physical value: F* = 186 Mev . The flavour symmetry 
breaking (FSB) in the Lagrangian will be considered below. 

The Wess-Zumino term, which should to be added to the action, is given by 



where ft is a 5-dimensional region with the 4-dimensional space-time as its boundary and 
where is a 5-dimensional analogue of = U^d^U. As is well known, this extra term does 
not contribute to the static masses of classical configurations, but it defines important 
topological properties of skyrmions || |[J and plays an important role in their quantization 



where r, 6 and ip are polar coordinates, /(r) is a radial profile function which has to be 
determined numerically and P(0, ip) is a projector involved in the harmonic map ansatz. As 
shown in |3| this projector is given by 



where F is a 3-component vector, whose entries are polynomials in z = tan(^)e lv . The 
largest degree of the polynomial gives the baryon number B of the final 577(3) Skyrme field 
configuration. For an SU(2) embedding the approach is similar except that this time we 
put 



where U 2 is an U{2) matrix determined in an analogous way as (3) except that the vector F 
has only two components. 

The fields F in Eq. (4) are so chosen that the configuration (3) has the smallest 
energy; this provides us with the harmonic map approximation to the real minimal energy 
static solution of the Skyrme model. We have then taken the expressions for F determined 
in and used the corresponding U in (3) as an initial field in our minimization program. 
For B = 2 F= (l,V2z,z 2 ) T , for B = 3 F = (l/y/2, 1.576z, z 3 ) T , for B = 4 F = (1, 2.72z 2 , z 4 ) T , 
for B = 5 F = (4.5z 2 , 2z 4 + 1, z 5 - 2.1 z) T and for B = 6 F = (fcz 3 , 1 - 3z 5 , z 6 + 3z f with k = 7.06 @ . 

We have performed many minimizations of the field configurations involving baryon 
numbers 3,4,5 and 6. In all cases the initial 577(3) harmonic map fields had energies some- 
what higher than the corresponding harmonic map embeddings. 

Our 577(3) minimization program, which uses the parametrization of the SU(3) field 
in terms of two, mutually orthogonal, complex unit vectorsj] was then used to minimize the 
energy further. The constraints of the orthogonality and of the vectors being of unit length 
were replaced by the introduction of extra positive contributions to the energy (with large 
coefficients) - the so called "penalty terms" . 

The program itself minimized the energy using a mixture of a finite element and 
finite step methods, and we varied the coefficients of the penalty terms so as not to be 
trapped in a local minimum. In each case the final energy of the field configuration was 

lr The authors thank W.K. Baskerville for making them aware of this parametrisation 
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lower by a few hundreds of MeV from the energy of the harmonic approximant. However, 
due to the smallness of the lattice the baryon number was also somewhat lower; so our 
values in the Table were obtained by a linear extrapolation: M = M + AM, where AM = 
a* (B - B obs ) and where B is the baryon number and B obs is its value on the lattice, a was 
determined by looking at various minimizations (with different values of the coefficients of 
the penalty terms). 

To check the stability of the program we also performed a series of minimizations 
when the initial field was a mixture of the SU(3) harmonic map field and of the embedding. 
When the mixture was close to the embedding or to the SU(3) harmonic map the minimiza- 
tion program took it down to these fields. When the mixture was far from either of these 
special fields it evolved to a new configuration of higher energy than either of the special 
ones. Although some of these new configurations can be numerical artifacts it is clear that 
the spectrum of static solutions of the SU(3) model is very complicated, with most states 
having energies larger than the energies of the embeddings or the fields derived by the 
harmonic map ansatz. Hence our estimates in the Table have some physical justification. 

In the Table we present some of our results. It is difficult to assess their accuracy; 
it is probably within 1% for the masses and several % for other quantities. In all numerical 
minimizations we have worked on a small tetrahedral lattice and hence we had a small 
"leakage" of both the energy and of the baryon number. 
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Table. The values of the masses M c i in GeV, the mass term AEM.t. (m MeV), the strangeness content Cs 
and the moments of inertia (in GeV~ x ) for the SU(3) projector configurations. The quantities for B = 1 
hedgehog and the masses of SU(2) embeddings M^ 2 , in GeV are given for comparison. The components 
of the Wess-Zumino term W Z% and WZg, which are different from zero for the SU(3) projector ansatz are 
also shown. AE in the last column is the quantum correction (in MeV) due to zero modes for the lowest 
energy state. The parameters of the model are — 186 MeV, e = 4.12. 



As can be seen from this Table, the moments of inertia are pairwise equal, except 
for the case B = 5 when 4 ^ 9 5 . These equalities are a consequence of the symmetry 
properties of our multiskyrmion configurations. 

It should be noted that the energy of the B = 4 configuration is very close to that of 
the B = 4 toroidal soliton obtained in I988jll| within the axially symmetrical generalization 
of the ansatz by Balachandran et al. 
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Within the accuracy of our calculations we 
cannot be certain which configuration has the lowest energy in the chirally symmetrical 
limit. For 5 = 6 the torus-like SO(3) configuration has energy considerably higher than the 
energy of the SU(3) projector ansatz presented in the Table. 

Let us add here that the configuration of a double-torus form was considered also 



within the SU(2)- version of the model several years ago [14|. It has the energy somewhat 
higher than the energy of the known minimal energy configuration. 



3. Following [10, 11, 15] we consider the contribution of the Wess-Zumino (WZ) term 
which defines the quantum numbers of the system in the quantization procedure. Its ex- 
pression is given by (2) above. As usually, we introduce the time-dependent collective coor- 



dinates for the quantization of zero modes according to the relation: U(r,t) = A(t)Uo(r)A ! <(t). 
Next we perform some integration by parts and rewrite the expression for the WZ-term 
as: 

. J TrA^ A{R a Rf}R~ ( + L a L L 7 )d 3 x, (6) 
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where L a = U^d a U = iL k , a \ k and R a = d a U G ul = U L a U^, or 

f k j^u k WZ k d*x = k f^u k Lr, (7) 



- wz 



k=l fe=l 



with the angular velocities of rotation in the configuration space defined in the usual way: 
A^A = -^u k \ k . Summation over repeated indices is assumed here and below. The functions 
WZ k can be expressed through the chiral derivatives L k : 

WZi = WZ? + WZ? = (Ri k {U ) + 5 lk )WZ£, (8a) 

i, k = 1, ...8, and are given by [ fTo] 

WZ? = -{L U L 4 L 5 + L 6 L 7 ) - {L 2 L 3 L 8 )/VS - 2(L 8 , L 4 L 7 - L 5 L 6 )/V3 

WZ? = -(L 2 ,L 4 L 5 + L 6 L 7 ) - (L 3 LiL 8 )/V3 - 2(L 8 , L 4 L 6 + L 5 i 7 )/%/3 
T^Z 3 L = ~(L 3 .L 4 L 5 + L 6 L 7 ) - [L 1 L 2 L & )/Vz - 2(L 8 , L 4 L 5 - L 6 L 7 )/V3 
WZ? = -(L 4 , L 4 L 2 - L 6 L 7 ) - {L 3 L 5 L S )/V3 + 2(L 8) L1L7 + L 2 £ 6 )/%/3 
WZ 5 L = - (L 8s iiL 2 - ^£7) + (L 3 i4L 8 )/V3 - 2(L 8 , LiL 6 - L 2 L 7 )/Vs 
WZ? = (L 6 , LiL 2 + L 4 L 5 ) + (£ 3 i7i 8 )/V3 - 2(2 8 , L^b - £ 2 L 4 )/V3 
T^Z 7 L = (L 7 , LxL 2 + L 4 L 5 ) ~ (£ 3 i 6 £ 8 )/V3 + 2(L 8 , + L 2 L 5 )/V3 

WZ? = ->/3(LiL 2 i 3 ) + (i8i 4 i 5 ) + {L 8 L 6 L 7 ), (9) 
where (L 1 ,L 2 L 3 ) denotes the mixed product of vectors Li, L 2 , £3, «e (Li,L 2 L 3 ) = (L\ ■ L 2 A L 3 ) 
and L 3 = (L 3 + V3L s )/2, L 8 = (V3L 3 - L 8 )/2, L 3 = {-L 3 + V3L 8 )/2, L 8 = (V3£ 3 + L 8 )/2 are the 
third and eighth components of the chiral derivatives in the (u, s) and (d, s) S'[/(2)-subgroups. 
Rik(U ) = \TrXjJJ Q \ k ul is a real orthogonal matrix, and WZ? are defined by the expressions 
(9) with the substitution L k — > R k . Relations similar to (9) can be obtained for WZ 3 and 
WZ 8 ; they are analogs of WZ 3 and WZ S for the (u, s) or (d, s) SU(2) subgroups, thus clarifying 
the symmetry of the WZ-term in the different 577(2) subgroups of 577(3). 

The baryon number of the 577(3) skyrmions can be written also in terms of Li in a 
form where its symmetry in the different 577(2) subgroups of 577(3) is more explicit: 

B = ~i / (( L l> L 2L 3 ) + (L 4 ,L 5 £3) + (£6,£7^^ (10) 

The contributions of the three SU(2) subgroups enter the baryon number on an equal 
footing. In addition, mixed terms corresponding to the contribution of the chiral fields 
from different subgroups are also present. 

The results of calculating the V^Z-term according to (9) depend on the orientation 
of the soliton in the S77(3) configuration space. The Guadagnini's quantization condition 
JT0] was generalized in [|T5] for configurations of the "molecular" type to 



2 dL wz _ N c B(l-3C s ) 
71 du 8 



Y ™ m = ~ I " c±Jy \ J ^ b ? (11) 



where Yr is the so-called right hypercharge characterizing the SU(3) irrep under consider- 
ation, and the scalar strangeness content Cs is defined in terms of the real parts of the 
diagonal matrix elements of the matrix U: 



C S 



< 1 - ReU 33 > 



< 3 - Re(U n + U 22 + U 33 ) >' 



(12) 



where <> denotes the averaging or integration over the whole 3-dimensional space. When 
solitons are located in the (u,d) SU{2) subgroup of SU{2>) only L x , L 2 and L 3 are different 
from zero, Cs = 0, WZ^ and WZ$ are both proportional to the B-number density, and the 
well known quantization condition JlOj takes place 



Y R = N c B/3. 



(13) 



The interpolation (11) does not work so well for configurations we consider here. 

The expression for the rotation energy density of the system depending on the 
angular velocities of rotations in the £17(3) collective coordinate space defined in Section 2 



can be written in the following compact form [15 



F 2 
32 



Lrot — ( &1 



16e 2 



(Sl2+S 45 ) +(s 4 5+s 67 ) +(s 67 -si 2 ) + - (2§ 



S46 



s 57 ) 2 



(2s 23 + Sir - s S6 ) + 



-(2s 34 + §16 - S 27 ) 2 + (2s 35 + §17 + S 26 ) 2 + (2§36 + §14 + S*25) 2 + (2§37 + §15 " S 2 ±Y 



(14a) 



or 



L ro t = Vik^iUJk/2 = Vikgug m kUiui m /2 



(146) 



Here s lk = u>iLk ~ u k Li, i,k = 1,2. ..8 are the 51/(3) indices, and s 34 = (§34 + \fis m )/2, s 35 = 
(s 35 + V3s 8 5)/2, s 36 = (-§36 + V3s 86 )/2, s 37 = (-§37 + v / 3§*87)/2, ie similarly to the definitions of 
L 3 and L$. gu are given in (15). To get (14) we have used the identity: s ab s cd - s ad s cb = s ac s bd - 
The formula (14) possesses remarkable symmetry relative to the different SU{2) subgroups 
of SU(3). The functions L$ or L 8 do not enter (14) nor the expression (10) for the baryon 
number density. The functions cDj are connected to the body-fixed angular velocities of 
SU(3) rotations by the transformations u = UqU>U - w, or, equivalently 



ft* = (Rik(U^) - 8 ik )u) k = g k iU) k . 



(15) 



Here R ik (V^) = R ki {V) is a real orthogonal matrix, i,k=l, ...8, and u>f = 2(wf - R k i(U )Lu k uji). 

The expression for the static energy can be obtained from (14) by means of the 
substitution S>i — * 2Li and s lk ~> 2n ik , |15| with n ik being the cross product of an L k ie 
il lk = Li A L k - From (14) we have then the inequality |15[ 



E stat - M.t. > 37r 2 B(J^/e), 



(16) 



which was obtained first by Skyrme for the SU(2) model. 

Eight diagonal moments of inertia and 28 off-diagonal ones define the rotation en- 
ergy, a quadratic form in ujiuj k as follows from (14) and (15). The analytical expressions for 
the moments of inertia are too lengthy to be reproduced here. Fortunately, it is possible 
to perform calculations without explicit analytical formulas, by substituting (15) into (14). 



For configurations generated by SU(3) projectors the Lagrangian of the system can 
be written in terms of angular velocities of rotation and moments of inertia in the form 
(in the body-fixed system): 

L rot = ^1(ojI+ojI) + + ^-coj + ^uj 2 5 + ^(cu 2 6 +u 2 7 ) + ^col + e 38 co 3 u 8 + W Z 3 U 3 + WZ 8 U 8 . (17) 

After the standard quantization procedure the Hamiltonian of the system, 
H = u>idL/duji - L, is a bilinear function of the generators Ri of 5*7(3) rotations: 

= R 2 + R 2 (R 3 - WZ 3 ) 2 _Rf_ Rj + R 2 

20i 8 2D 38 26 4 26 5 26 6 

+e 3 {Rs ~^ Zs)2 -^(Rs- WZ 3 )(R S - WZ 8 ), (18) 

ZL> 3 8 L> 38 

where D 38 = 636s - ©is- For the states belonging to a definite 51/(3) irrep the rotation 
energy can be written in terms of the second order Casimir operators of the SU(2) and 
5*7(3) groups: 

_ N(N+1)-R 2 U(U + l)-Rj ua V(V + l)-Rj dS i 
rot ~ 26i + 26 4 29 6 + 

s 8 (r 3 -wz 3 ) 2 q 3 (r 8 -wz 8 ) 2 e 38 (iff--Rg)(e B -8 4 ) 

+ 2^ + ^ 7^ (i?3 - - ^ + 46^ (19) 

with e 4 = 2e 4 e 5 /(e 4 + e 5 ). 

R3,us = R3/2 + VSR 8 /2, R 8 , us = V3R 3 /2 - R s /2, 

Rs,ds = -^3/2 + V3i?8/2, Rs.ds = V3i? 3 /2 + iZs/2, (20) 

with R 3 = R 3 . ud , R 8 = R 8 ,ud', -R 3 ,afc, Rs,ab denote the 3 rd and 8 th generators of the (a, 6) 5*7(2) 
subgroup of 5*7(3). 

TV, *7 and V are the values of the right isospin in the (u,d) 5*7(2) subgroup, and so 
called *7-spin and y-spin. They are connected with the second order Casimir operator 
of the 5*7(3) group C 2 (5*7 3 ) = \{p 2 + q 2 + pq) + p + q; p,q being the numbers of the upper 
and low indices in the tensor describing the 5*7(3) irrep (p,q) via the symmetric relation 
C 2 (5*7 3 ) = N(N + 1) + *7(*7 + 1) + V(V + 1) - (Rj + R%)/2. The hypercharge Y = 2R 8 /V3. The 
terms linear in the angular velocities present in the Lagrangian due to the Wess-Zumino 
term cancel in the Hamiltonian. 

For a state with the lowest energy which belongs to the 5*7(3) singlet with (p, q) = 
the quantum correction simplifies to: 

ae = -l—(e s wz 2 + e 3 wz 2 - 2e 38 wz 3 wz 8 ). (21a) 

2L) 38 

The quantities Z? 38 , 6 3 + 6 8 , WZ 2 + WZ 8 do not depend on the choice of the particular 
subgroup of 5*7(3) group, therefore Eq. (21) also is invariant, as one expects from general 
arguments, since 

Q & WZl + Q 3 WZl - 20 38 WZ 3 WZ 8 = (6 8 + 3 ){WZ 2 + WZ%) - <d 3 WZ 2 z - <d 8 WZ\ - 2Q 38 WZ 3 WZ 8 . 
The generalization of (21a) to an arbitrary case is 

AE min = QT.^WZiWZj/2, (21b) 



where O" 1 is the matrix inverse to the tensor of inertia 9 y . Eq. (216) is valid for all cases 
except the degenerate case when DetQij = 0. Note, that the case of the SU(2) embedding just 
corresponds to the degenerate case, with 9 8 = 0, and this leads to the rigorous quantization 
condition [jTOj . 

As can be seen from the Table, for the configurations we consider here and for our 
choice of the subgroups 9| 8 < 6 3 9 8 , so a small correction due to 9 38 can be neglected, and 
AE ~ WZ$/(2<d 3 ) + WZ$/ (20$). Since the moments of inertia are of the order of magnitude 
— lOGeV^ 1 and both WZ 3 and WZ 8 are not greater — 1, the quantum correction for the 
lowest states does not exceed several tens of MeV: see the Table. 

The quantum correction from the SU(3) zero modes which should be added to the 
classical energy of the soliton located originally in the (u,d) SU{2) subgroup equals to 
AE = 35/49 F - 365 MeV [TIE 



Note, that for the 50(3) solitons considered in O, O the WZ-term is equal to zero 



and for this reason the lowest quantized state was an SU(3) singlet without any quantum 
correction to its mass. 

4. The FSB part of the mass terms in the Lagrangian density which defines, in 
particular, the mass splittings inside 5C/(3) multiplets is defined, as usual, by 

L FSB = F D m D- F >l Tr{l _ V3A 8 )(C7 + C/t _ 2 ) + F ^~^ Tr{\ - s/3\ s )(Uy» + l^rf) (22) 

Here m D is the mass of K, D or B meson. The ratios F D /F^ are known to be 1.22 
and 1.7 ±0.2 for, respectively, kaons and D-mesons. The L F sb given by (22) is sufficient 
to describe the mass splittings of the octet and decuplet of baryons within the collective 
coordinate quantization approach 

The contribution to the classical mass of the solitons from the 5t/(3) projector ansatz 
which comes from FSB part of the Lagrangian (22), considered as a perturbation, is not 
small even for strangeness: 

AM = 2C s AE M . t [^^-l^ (23) 

where AE M .t is the FS mass term contribution shown in the Table. Numerically, it equals 
1.13 GeV for B — 2 and 3.2 GeV for B = 6. Thus it makes sense to include the FSB part of 
the mass term into the minimization procedure. 

We do not consider here the quantum corrections due to the rotations in the ordinary 
space because they are explicitily zero for the lowest states with even B ( J = o) , or small 

Our investigations have shown that the space of local minima for the 5f(3) skyrmions 
has a rather complicated form. In addition to the known local minima - the SU(2) em- 
beddings, 50(3) solitons and skyrmion molecules - there are also other local minima best 
approximated by the 5C(3) projector ansaetze [|3], [|. 

It is possible to quantize these configuratons and to estimate the spectrum of states 
by means of the collective coordinate quantization procedure. The total sum of the classical 
mass and of the zero modes quantum corrections can be the smallest for some baryon num- 
bers. However, to draw the final conclusions one needs to calculate the Casimir energies of 
solitons which are essentially different for solitons of different form and size. This is a very 



complicated problem solved approximately [19| -[21] only for the B = 1 hedgehog configu 



rations. For the case of SU{2) embeddings it was, however, possible to draw conclusions of 
physical relevance under the natural assumption that the unknown Casimir energy, or loop 
corrections, cancel in the difference of energies of states with different flavours, and that 



the states with u, d flavours can be identified with ordinary nuclei [^2] . But this assumption 
may be incorrect for the states generated by the SU(3) projector ansatz. 

We are thankful to W.K. Baskerville, T. Ioannidou, B. Piette and P.M. Sutcliffe for 
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